Spatial structure of a vortex in low density neutron matter 
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We study in a fully selfconsistent approach the structure 
of a vortex in low density superfluid neutron matter. We 
determine that the matter density profile of a vortex shows a 
significant depletion in the region of the core, a feature never 
reported for a vortex state in a Fermi superfiuid. 

PACS numbers: 26.60. +c, 21.65.+f, 97.60. Jd 



Almost thirty years ago Anderson and Itoh [1] put 
forward the idea that vortices should appear in neutron 
stars and that they can also get pinned to the solid crust. 
They argued that the "star-quakes," observable on Earth 
as pulsar "glitches," apparently the only evidence so far 
that solid matter exits in the universe apart from plan- 
ets, are caused by the vortex de-pinning. This idea and 
its various implications have been examined by numerous 
authors, see Refs. [2] and further references therein, but a 
general consensus does not seem to have emerged so far. 
In spite of this, there does not seem to be a doubt in any- 
body's mind that a significant part of the matter in a neu- 
tron star is a superfluid of one kind or another. In partic- 
ular, for densities less than nuclear saturation densities 
one expects neutrons to form ^5*0 Cooper pairs, similar to 
ordinary electrons in superconductors. Even though at 
low neutron densities the kinetic energy dominates over 
interaction energy, since this is mainly attractive. Cooper 
instability occurs. Neutrons become a Fermi superfluid 
and vortices can appear. In a low density neutron mat- 
ter various nuclear-like objects [3] exist as well, a state 
of matter often referred to as "pasta," "meat balls" and 
"Swiss cheese" phases. The nuclei-like objects in some of 
these phases are significantly denser than the surround- 
ing low density neutron matter and positively charged 
and are expected to form a rather stiff Coulomb lattice. 
However, other alternatives (liquid crystals, amorphous, 
disordered or heterogeneous state) have been discussed 
in the literature as well [4]. In any case, these inho- 
mogeneities are the "islands" on which vortices can get 
"pinned down," at least temporarily. The estimates of 
the de-pinning energy available in literature are rather 
crude. The profile of a vortex in neutron matter is typ- 
ically determined using a Ginzburg-Landau equation, 
which is expected to give mostly a qualitative picture 
and its accuracy is difficult to estimate. Surprisingly, 
there exists only one microscopic calculation of a vortex 
in low density neutron matter [5]. In regular supercon- 
ductors the pairing energy is a relatively small quantity, 
when compared to the Fermi energy, and the presence or 



absence of the pairing field is not expected to lead to any 
noticeable changes in the electron density [6]. For this 
reason is typical in the case of electrons to neglect the 
mean-field (more exactly, to assume it a constant), and 
only account for the presence of the magnetic field in- 
duced by the electron super-flow. It is known that inside 
a vortex core one can have bound states within the pair- 
ing gap [7], basically due to the phenomenon known as 
Andreev reflection. The presence of these localized states 
inside the vortex core leads to a number of observables, 
in particular their presence modifies the specific heat at 
low temperatures. This aspect could prove crucial for the 
thermal emission from neutron stars, but apparently so 
far this aspect has not been considered. A small change 
in the profile of the matter density and of the mean-field 
was recently reported in dilute atomic Fermi gases, but 
it was not considered noteworthy [8], see also Ref. [5]. 

We report here on a fully microscopic analysis of a 
vortex in pure neutron matter and but we do not con- 
sider explicitly so far an impurity, on which such a vor- 
tex might get pinned. The ratio of the pairing gap to 
the Fermi energy in low density neutron matter is larger 
than in other systems and as we have found, the matter 
distribution is influenced appreciably by the presence of 
a vortex. It is natural to expect that a major change 
in the density profile of a vortex is bound to affect the 
magnitude and the character of the pinning mechanism. 
This particular aspect, apparently never reported in lit- 
erature until now, can have significant consequences on 
the physics of vortices in neutron stars. 
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FIG. 1. The energy per particle as a function of neutron 
density for the non-interacting neutron gas, interacting with 
and without pairing correlations taken into account, for the 
case of homogeneous matter distribution. 

In order to describe selfconsistently neutron matter. 
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FIG. 2. The neutron density and neutron pairing gap as 
a function of the neutron chemical potential in the case of 
homogeneous neutron distribution. 



we use a nuclear energy density functional approach 
(NEDF), which works surprisingly well in finite nuclei 
[9-11]. This NEDF has as the main part the contribu- 
tion describing normal nuclear matter [12], which was fit- 
ted by Fayans to reproduce the microscopic many body 
calculations of Refs. [13] for the case of homogeneous 
matter distribution. Fayans has supplemented this in- 
put with a phcnomcnological gradient correction energy 
density (which we take into account here) as well as with 
a spin-orbit energy density (neglected here). A separate 
contribution to the energy density functional should de- 
scribe pairing correlations. It is well known that even 
in very dilute systems, the induced interactions play a 
significant role in determining the magnitude of the pair- 
ing gap, drastically reducing it by a factor close to 2.2 
when compared to a naive BCS calculation with bare 
interaction [14-16]. Neutron matter becomes dilute in 
the sense that p\a\^ < 1 (where a is the nn-scattering 
length) only at densities smaller than « 6 x 10~^ fm~^. 
The range of neutron densities we are considering here 
are significantly larger and the arguments presented in 
Refs. [14-16] strictly do not apply. Over the years a 
number of investigators arrived at qualitatively similar 
conclusions, even though the specific corrections consid- 
ered often varied from one set of authors to another [17] 
and the neutron matter was certainly not dilute in the 
sense mentioned above. One should note also, that in nu- 
clei [18] the effect seems to be the opposite, a conclusion, 
which for the wrong reasons, agrees with the analysis of 
Ref. [15] for the case when the number of fermion species 
is four (two for spin and two for isospin). All analyses 
however seem to agree that the pairing gap in neutron 
matter as a function of the Fermi momentum has a max- 
imum of approximately 1 MeV for a Fermi momentum of 
kp ~ 0.8 fm~^. and that the gap vanishes for Fermi mo- 
menta larger than ft! 1.5 fm~^. This shape of the pairing 
gap could be reproduced rather satisfactorily with the 
formula [19] 
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where fc^T- is the Fermi wave vector, m the nucleon mass 
and S{kp) is the ^5*0 nn-phase shift. Following the argu- 
ments given in Refs. [14,15], we have changed the corre- 
sponding formula for the gap given in Ref. [19] by includ- 
ing the effect of induced interactions as an additional fac- 
tor of l/(4e)i/3 ft 0.4 5. Since for the neutron densities we 
are considering here, the effective range approximation to 
the Nijmegen ^5*0 nn -phase shift is essentially exact [20], 
we have parametrized the pairing gap in homogeneous 
neutron matter using this approximation. The resulting 
pairing gap is well within the theoretical uncertainties, in- 
cluding various corrections as discussed in Refs. [17] and 
references therein. For each value of the neutron den- 
sity we have determined a corresponding value for the 
"bare coupling constant" g{p) as discussed in Refs. [10]. 
This completes the construction of the NEDF, in com- 
plete analogy with the construction of the normal part 
of the electron EDF [21], see Figs. 1-2 for details. 
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FIG. 3. The neutron density as a function of the distance 
from the vortex axis r. Some of these density distributions 
have boon rescaled by the factor indicated above the corre- 
sponding curve. 

The rest of the technical details were described pre- 
viously in Refs. [9-11]. We note only that the normal 
and anomalous densities were evaluated using the tech- 
nique initially described in Ref. [22], and consequently 
all the boundary conditions were taken into account ex- 
actly and there was no need to either enclose the vortex 
in a cylinder or to introduce periodic boundary condi- 
tions along the vortex axis. Moreover, there is no need 
to diagonalize any matrices and the fact that the single— 
quasiparticle spectrum has (or has not) a number of dis- 
crete eigenvalues does not require any change in the nu- 
merical algorithm. Since we consider a straight vortex, 
the quasiparticle wave functions in cylindrical coordi- 
nates r = (r, <j), z) have the simple structure 
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Ua{r) exp [in4) -\- ikz] 
Va{r) exp [i{n — + ikz] 
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FIG. 4. The mean-field neutron potential as a function of 
the distance from the vortex axis. Some of these potentials 
have been rescaled by the factor indicated above the corre- 
sponding curve. 

where a labels the quasiparticle states, n is integer and 
k is the wave vector of the quasiparticle state along the 
vortex Oz-axis and a = (a, k, n). (Note that wc use r for 
the 3d-coordinate and r for the radial cylindrical coordi- 
nate.) The mean-field depends only on radial coordinate 
r (measured from the vortex axis), while the pairing field 
has the structure A(r) exp(z(/)). 

Our results are summarized in Figs. 3-6 where we 
show the vortex density, mean field, pairing field and 
velocity profiles, as functions of the distance from the 
vortex axis, for several values of the neutron chemical 
potential. The most unexpected feature of these profiles 
is the prominent depression of the matter density in the 
region of the vortex axis. In the case of superfluid dilute 
bosons at the vortex axis the density practically vanishes. 
Except for the barely visible features in the density profile 
of a vortex reported in Refs. [5,8] but not commented on, 
such features have not been discussed in literature in the 
case of a Fermi system. The asymptotic density changes 
roughly by a factor of forty from the lowest to the highest 
values of the chemical potential plotted in these figures 
and the Fermi momentum changes by approximately a 
factor of three. The size of the inhomogeneity in the 
matter distribution and the mean-field U{r) are governed 
by the asymptotic value of the Fermi wavelength, while 
the spatial profile of the pairing field A(r) is controlled 
mostly by the coherence length, ^ = ep/^kp ^ 1/kp 
(where kp, A and ep are the asymptotic values for the 
Fermi wave vector, pairing gap and Fermi energy). The 
density gradients are noticeable and the gradient correc- 
tion terms in the NEDF, even though they are small, are 
non-negligible. The density depletion drives the mean- 
field to become less attractive in the vortex core, which 
in its turn, due to the self-consistency, stabilizes this 
structure. The density depletion is likely so significant 
in neutron matter because the pairing field is relatively 
stronger than in the case of electrons. The magnitude 
of the pairing correlations is characterized by the ratio 



A{r)/ep{r) = A(r)/(/i — I7(r)). In the vortex core region 
the behavior A(r) — > as r ^ can be partially com- 
pensated by an increasing U{r), which induces in its turn 
a density depletion. The density and pairing field pro- 
files (and to a lesser extent the mean-field as well) show 
some rather faint Friedel-like oscillations, which can be 
attributed to the presence of discrete states inside the 
vortex core [5-8] . 
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FIG. 5. The pairing gap as a function of the distance from 
the vortex axis. 

The presence of the vortex implies a net fiow of the 
neutron superfluid around the vortex axis with the ve- 
locity given by 

ill d 

^^^^^^^ = -d;^^'"^'"^'''^'"^'"^' 

where = (— y, x,0)/r, the summation over a should be 

interpreted as a sum or integral when appropriate, over 
all quasiparticle states with energies Ea > and p{r) is 
the neutron number density. Naturally Vs(0) = 0. If all 
neutron matter would be involved in such a super-flow 
around the vortex axis, the superfluid velocity would be 
given by 

Vv(r)e^=^(-y,x,0), (4) 

which corresponds to Osanger's quantization condition 
[23] per Cooper pair 

i-£v.(r)e,.dr = A, (5) 

and where the contour C is any closed curve, around 
the vortex axis. While at relatively large distances from 
the vortex axis the neutron super-flow indeed approaches 
this "classical" limit, that is not the case at distances 
smaller or of the order of the coherence length ^. The 
interpretation of this behavior is relatively simple. Since 
the pairing field vanishes at the vortex axis, the neutron 
matter in the immediate neighborhood is not fully super- 
fiuid. Only the fraction fs{r) = Vs(r)/Vv(r) of the neu- 
trons are entrained into the superfluid flow. This fraction 
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vanishes at the vortex axis and approaches unity far from 
the vortex core, essentially in a monotonic fashion. Vs(r) 
is always significantly smaller than the Fermi velocity. 
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FIG. 6. The neutron current velocity Vs{r) (in units of 
speed of light) as a function of the distance from the vor- 
tex axis. The thick dashed line is a plot of the "classical" 
current velocity Vv{r) = h/2mr. 

The existence of a strong density depletion in the vor- 
tex core is going to affect appreciably the energetics 
of a neutron star crust. One can obtain a gross esti- 
mate of the pinning energy of a vortex on a nucleus as 
Epin = [£{Pout)Pout - e{pin)pin]V, whcic e{p) is the en- 
ergy per particle at density p, see Fig. 1, and pout 
are the densities inside and outside the vortex core and 
V is the volume of the nucleus. Naturally, this simple 
formula does not take into account a number of factors, 
in particular surface effects and the changes in the ve- 
locity prohle and the pairing held. These last contribu- 
tions were accounted for (with some variations) in the 
past [1,2]. However, if the density inside the vortex core 
and outside differ significantly one expects to be the 
dominant contribution. In the low density region, where 
£{Pout) Pout /siPin) pin IS largest, see Fig. 4, one expects 
a particularly large anti pinning effect (Ep^^ > 0). The 
energy per unit length of a simple vortex is expected to 
be significantly lowered when compared with previous es- 
timates [1,2] by w [e{pout)Poui - £{pin)Pin]'n:R^ ■, where R 
is an approximate core radius, see Fig. 3. How a sig- 
nificant density depletion found here affects the neutron 
star properties in the crust will be addressed in our future 
work. 
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